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I. INTRODUCTION 



One of the most interesting features of QCD (quantum chromodynamics) is the transition from the confinement to 
the deconfinement phase. In this new state of QCD, quarks and gluons confined in tlic hadron at zero temperature 
move freely when the system reaches a sufficiently high temperature. The quark gluon plasma (QGP) was realized 
at high temperature in the early universe, and is expected to be produced in heavy-ion collision experiments at the 
CERN Super Proton Synchrotron (SPS), BNL Relativistic Heavy Ion Collider (RHIC), and CERN Large Hadron 
Collider (LHC). Thus it is an urgent task to accumulate theoretical knowledge about the QGP. 

The massless gluon in the QGP medium is changed into a dressed massive gluon after quantum corrections. The 
screening effect is characterized by a mass pole of the propagator and is closely related to thermal QCD phenomenology. 
One example is a screened heavy quark potential, which is frequently discussed in relation to J/V' or T suppression. 
For calculations of jet quenching, which might be a fingerprint of a QGP, a model including the electric and magnetic 
masses has been proposed . A nonperturbative quantitative study in the vicinity of Tc and up to several times 
is of great importance for understanding QGP physics. 

The thermal field theory is the most basic method of studying the QGP and has provided many informative 

observations. At zero temperature, many calculations based on pcrturbative QCD have described experiment results 
[3 and there is no doubt that QCD is a theory of strong interaction. It is natural to employ the pcrturbative approach 
to thermal QCD. Because of the asymptotic freedom at high temperature, the coupling constant is expected to become 
small enough to carry out the perturbations. In such a high energy state, quarks and gluons must behave as an ideal 
gas; yet this simple consideration is spoiled by an infrared divergence 6], which is known to bear a hierarchy on 
the energy scale in the QGP system [HQ. One usually defines l/T as a perturbative length scale, and furthermore 
1/gT must be introduced as an electric (Debye) scale, whose influence appears as a Yukawa-type potential rather 
than a Coulomb-like one, and as a magnetic scale. In QCD, the magnetic mass, which cannot be accessed 

perturbatively, acts as a cutoff factor in the infrared problem and consequently becomes an essential element of 
thermal QCD. 

The QCD coupling constant strength near the critical temperature is still of the order of 1. Therefore perturbation 
theory is not applicable. However, these regions are currently being investigated with great interest in much theoretical 
and experimental research. After hard-thermal-loop (HTL) resummation was consistently formulated by Braaten and 
Pisarski, there were several improvements 0, for unsolved problems. Comparison of this method with lattice 
numerical data has been reported: one-loop HTL calculations of the free energy of a QGP are in good agreement with 
the lattice numerical result However, a recent two-loop HTL calculation indicates that it does not yet have 

adequate convergence As another approach, 3D reduction theory has also been widely studied and has yielded 
some promising arguments }T3|, but this method, which is defined only for the high temperature limit, cannot be 
applied to confinement/deconfinement physics. 

For reliable phenomenological analyses of high energy heavy-ion collisions, it is important to obtain information on 
the magnetic and electric masses of gluons nonperturbatively; a numerical study of lattice QCD as a first-principles 
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calculation should play an important role here. 

There have been many lattice studies of finite-temperature QCD, but only a few calculations of the electric and 
magnetic screening masses can be found in the literature other than for the case of color 5'C/(2) [11 d. The 
electric mass has been estimated from the Polyakov loop correlation functions to obtain the screened heavy-quark- 
antiquark potential jl^ J/Tfe .ISJ ■ 

The main aim of our study is to obtain reliable SU{3) electric and magnetic masses through large-scale lattice 
QCD simulations and to reveal their temperature dependence [T9l |. We also compare our numerical data with the 
predictions of leading order perturbation (LOP), HTL resummation, and other analyses. 

Since our mass extraction is based on measuring a gauge dependent gluon propagator, it is important to check the 
gauge invariance of both screening masses . This test is essential particularly for the magnetic mass because it 
has a poor definition in the frame of the perturbation. 

Gluons are essential ingredients in QCD dynamics, and QCD undergoes a phase transition from the confinement to 
the deconfinement phase when the temperature increases. Therefore, we expect gluon propagators to show different 
behavior in each phase, and their study provides information on confinement/deconfinement dynamics 21] . 

We measure gluon propagators, which depend on the gauge used, and therefore a gauge-fixing procedure is indis- 
pensable. However, gauge fixing on the lattice is difficult practically and conceptually. Usually, gauge fixing is carried 
out by the iterative technique |22| . and it is very time consuming. The conceptual difficulty is that the gauge is not 
uniquely fixed; this is known as the Gribov copy problem 23] . In order to overcome these difficulties, we adopt here 
a stochastic quantization with Zwanziger's gauge-fixing term ^24, 25, 26], instead of the path integral method. 

In this paper, we extend our previous results from the 20^ x 32 x 6 lattice simulations [l^ to the level that the 
present computer power can reach; we add more detailed values for the electric and magnetic masses and results 
at higher temperature from the larger lattice 32^ x 48 x 6 simulation. We also give a detailed description of the 
algorithm employed in this study. In Sec. II, we describe the stochastic gauge quantization together with the Gribov 
copy problem. The definitions of the gluon propagator and electric and magnetic masses are given. A large part of 
Sec. HI is devoted to our simulation results on the small lattice 20^ x 32 x 6. First we describe all input parameters 
of the simulation and the statistics needed to measure reliable gluon propagators. Then we show the gluon behavior 
and the electric and magnetic masses extracted from it. The gauge dependence check and temperature dependence 
for both screening masses are also given. Finally, we compare the numerical data with the perturbative argument, 
add the higher temperature result, and comment on the finite volume effect. Section IV gives the conclusions. 

The main part of the calculation was carried out on the SX-5(NEC) vector-parallel computer of RCNP(CMC), 
Osaka University. We used a parallel queue with 4, 8 and 16 CPUs and required about six months to complete this 
work. 



II. STOCHASTIC GAUGE FIXING ON LATTICES AND GLUON PROPAGATORS 



A. Lattice gauge action 

The lattice regularization scheme of QCD is the gauge invariant Euclidean theory which enables us to perform a 
nonperturbative calculation based on the Monte Carlo numerical technique . The simplest standard Wilson gauge 
action of lattice QCD can be defined from the continuum QCD as 

Sg = - l^"" Tr [U,{x)U^{x + v)UUx + /i)C/t(a;)] ), /3 = (1) 

■J 9 

Here a link variable, U^{x), stands for the SU{'i) color gauge field: 

[/^(a;) = e^s''"^''^"^), (2) 

where is the lattice spacing, i.e, the lattice cutoff, and represents the gauge potential of the gluon. In this 
study, we adopt the quenched lattice simulation (pure gauge QCD) without a dynamical quark effect, using Eq. 

B. Gauge fixing and Gribov copy on lattices 

We are interested in a direct calculation of the gluon propagator, and the extraction of electric and magnetic 
screening masses from it. Therefore we must fix the gauge of the gluon fields on the lattice, where the gauge 
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transformation is given by 



Uf,{x) uJ^x)U|_,{x)uJ{x + fi). 



LO stands for a gauge rotation matrix G SU{i) on the lattice. 

In this study, we focus on a Lorentz-type gauge, which is defined in the continuum as 



while in the discrete lattice theory 



A°(a;) EE ^ 2 Im Tr e{U^{x) - U^,{x - A)} = 0. 

Here t° is the SU{N) generator with the relation Tr[i"i''] = |<5"''- The above condition ISJ is equivalent to 

6J = 0, 

/ = E.,,, Re Tru;\x)U^{x)u:{x + A). 



(3) 



(4) 



(5) 



(6) 



Wilson j23| and Mandula and Ogilvie |2.3| suggested the following condition for the gauge fixing method on the lattice: 



Max.^/. 

The continuum version of Eq. Q was discussed in Refs. and [SC 



(7) 



Gauge fixing plane 

a A =0 




FIG. 1: Gribov [23 pointed out that a gauge is not uniquely fixed for non-Abelian theories, which is called the Gribov copy 
problem. This does not appear in the perturbative frame. The existence of the Gribov copy has been confirmed in several 
works [sllsllsl. 



The conditions © and Q are not equivalent to each other since there may be local maxima or minima of / which 
satisfy Eq. The gauge- fixing configuration cannot be uniquely fixed; this is called the "Gribov copy" ^23] and is 
illustrated in Fig. ^ If we study this problem using a numerical lattice simulation based on the iterative procedure, 
it is very difficult to find a true maximum. 



C. Stochastic gauge fixing 



Our approach to the gauge-fixing procedure is to use the stochastic gauge quantization instead of the Monte Carlo 
path integral. The stochastic quantization is based on the Langevin equation which introduces virtual time in addition 
to the Euclidean coordinate. Zwanziger introduced a gauge-fixing term as 

where Df{A) is a covariant derivative, t stands for Langevin time, and 77 is a Gaussian noise term. The second term 
on the right-hand side is a gauge-fixing term, a is a gauge parameter; a = corresponds to the Lorentz gauge and 
a = I to the Feynman gauge. 
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Monte Carlo Steps 




Langevin Steps 

FIG. 2: The figure illustrates the gauge-fixing procedure by the Wilson-Mandula-Ogilvie iterative method and by stochastic 
gauge fixing. Gray arrows indicate the gauge rotation for each algorithm. In the first algorithm, Monte Carlo update steps are 
performed without restriction on the d^iA^i — plane. When a gauge fixed configuration is needed, the gauge configuration 
is rotated to the gauge-fixed plane, df^A^ — 0, by iteration (2^. On the other hand, in the case of Langevin gauge fixing, 
configurations are updated by fluctuation around d^iA^ — 0. 



Mizutani and Nakamura [265 developed the lattice version of the stochastic gauge fixing. The link variables are 
rotated through the following gauge transformation depending on the virtual time: 



Ufj,{x, T + At) = Lj^x, r)e'-''^*°[/^(a;, t)ll;{x + fi, r). 



Here / stands for the force 



dS 



and the gauge rotation matrix is given by 



„i/3A"t°Ar/a 



(9) 



(10) 



(11) 



li uj = I, Eq. jnj is a lattice Langevin process. Gauge rotation, Eq. 0, with Eq. leads to the gauge fixing term 

as At — > 0. Equation means that the gauge rotation and Langevin step are executed alternately, as illustrated in 

Fig. [a 

In this stochastic quantization with Lorentz-type gauge fixing, A"^ fluctuates around the gauge-fixing plane A" — 0. 
For example, when we take a — 1.0 and At = 0.01, A° on 4^ x 8 with /3 — 6.0 behaves as shown in Fig. |21 We 
confirm that A"^ fluctuates arou " ' ' ' ved. 




800 1000 



FIG. 3: These data are produced on the lattice 4'' x 8 with Ar = 0.01(Langevin step interval), a = 1.0 (gauge parameter) at 
/3 — 6.0. The solid line means the average of ^{A"' (x)}'^ , while the dashed line is the value of A^(l). A" fluctuates around 
A" — and gauge configurations are sufficiently fixed on the dfj,Afi — plane. 
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There are two reasons for using the stochastic gauge-fixing method in this study. One is a practical issue. When 
we use the standard Wilson-Mandula gauge-fixing method (23], the iterative procedure is appUed to accepted gauge 
configurations for each Monte Carlo step in Fig. |21 Then the number of iterations is unpredictable, particularly for 
large lattices. On the other hand, in algorithm J^, we simultaneously repeat the step of update and gauge rotation in 
Fig. |2]and are free from the convergence problem of gauge fixing. Therefore we can estimate the CPU time precisely. 
Moreover, this algorithm, in which the gauge parameter a can be changed at will, is advantageous when testing gauge 
invariance. 

A conceptual problem is the Gribov ambiguity ^2^. The algorithm has a noteworthy feature, i.e., the second term 
of Eq. JSJ gives rise to a configuration such that 

fE.A2<oin n 

> out of n. ^ ' 



Here Vt stands for the Gribov region, 

Gribov region : VL = {A1{x)\d^,A1{x) = 0, -d^,Df > 0}. (13) 

That is, the stochastic gauge-fixing term is attractive (repulsive) inside (outside) the Gribov region 0,11^ if we start 
from the trivial configuration {yl^ = 0}. Although our algorithm may not completely eliminate copies, we conclude 
that it is a more effective method. 

Since the update algorithm described here is not as oooular as the Metropolis or pseudo-heat-bath method, we 
show the autocorrelation of the D that of the pseudo-heat-bath. 




FIG. 4; Autocorrelation of the Polyakov loops for the Langevin algorithm with the stochastic gauge fixing (Ar = 0.03) and 
the pseudo-heat-bath update algorithms as a function of the iteration steps. This calculation is done a,t 13 — 6.1 on the 8^ x 4 
lattice. 



D. Definition of gluon propagators 



We define the gauge field (x) in terms of the link variables as 

A^{x) = 2^^''ImTrr[/^(a;). (14) 

a 

We calculate massless gluon correlation functions with finite momentum, 

G^,^{px,Py,Pt,z) = {TTA^{px,py,pt,z)A^{-px,-py,-pt,0)). (15) 

In our study, to avoid a mixture of the longitudinal and transverse modes, we adopt the transverse conditions 
Pi ■ Ai = and measure the partially Fourier transformed propagator including one momentum p^ = 2t:/Nx: or 
Py — 2'KjNy We will obtain the gluon mass using a lattice energy momentum relation. This is different from 

other calculations [l^ , where the mass extraction from zero momentum propagator was studied. 
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TABLE I: We estimate a lattice cutoff and its temperature scale by using mainly QCD_TARO fit function |3: 
Tc ~ 256 [si. 



33. We adopt 



p 


a-'[GeV] T[MeV] T/T, 


P 


a-^[GeV] T[MeV] T /Tc 


5.8 


1.33 


222 


0.86 


6.4 


3.52 


586 


2.29 


5.90 


1.62 


270 


1.05 


6.5 


4.12 


690 


2.69 


5.95 


1.77 


295 


1.15 


6.6 


4.60 


767 


2.99 


6.0 


2.04 


340 


1.32 


6.7 


5.24 


874 


3.41 


6.05 


2.09 


349 


1.36 


6.8 


5.96 


994 


3.88 


6.1 


2.27 


378 


1.47 


6.9 


6.76 


1128 


4.40 


6.2 


2.64 


447 


1.74 


7.0 


7.64 


1274 


4.97 


6.3 


3.05 


509 


1.99 


7.1 


8.61 


1436 


5.61 



In thermal perturbative QCD the electric mass is defined from the temporal part of the gluon polarization 

tensor II^i^, while its spatial part is considered as the magnetic mass. Thus we can construct an electric propagator 
using a temporal one: 



2ir 27r 
G„(— ,0,0,z) + G„(0,— ,0,z) 



In the same way, a magnetic propagator is defined by spatial components: 



G..(0,^,0,z) + G,,(^,0,0,z) 



(16) 



(17) 



Because the screening mass is given as a pole of the denominator of the momentum space propagator, l/(p^ + m^), 
Ge{z) and G„i{z) are expected to behave as exponential damping functions in the z direction on distances, Z ^ 1/T 
with the masses, 



Ge{m) (2) 



)(P)2 



(18) 



III. RESULTS 



A. Simulation parameters 



We use mainly the lattice of the size N^NyN^Nt — 20^ x 32 x 6 which satisfies the condition ^ 3A^i. On 
this lattice the long range area corresponds io z ^ 1/T and thus reliable information for screening physics may be 
produced. 

We summarize lattice cutoff values and the corresponding temperature in Table Varying /3, i.e., the lattice 
coupling constant, we change the temperature T = 1/Nta. The pure gauge lattice with Nt — 6 has critical /3c ~ 5.89 
and we adopt Tc ^ 256 36] as the critical temperature. To estimate the lattice cutoff value, we employ the results in 

Refs. [Mllg. 

In addition, we prepare a lattice of the size 32^ x 48 x 6 in order to investigate the finite size effect of screening 
masses, and to obtain them at higher temperature. 



B. Necessary statistics to obtain reliable gluon propagators 



We observe a large fluctuation of the gauge propagators, particularly at long distances, i.e., we suffer from a long 
autocorrelation time. We show in Fig. |Slthe typical behavior of gluon propagators G{z) as a function of the Langevin 



step. In order to analyze the gluon propagator and measure the screening masses we require 0.2 
the typical number of simulation data. 



0.4 X 10^ steps as 
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T/T =2.69, Ai:=0.05, a=1.0 
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FIG. 5: Typical gluon propagator behavior as a function of Langevin steps on the lattice 20 x 32 x 6. We find that fluctuations 
of Gxx{12) are much larger than those of Gxx{6)- In order to investigate the screening effect, the long range contribution should 
be adopted. Consequently, we need a large number of statistics. 



C. Gluon propagators 

Although the gluon propagator itself is gauge dependent, it gives us some insights into the gluon dynamics of the 
confinement /deconfinement physics . 

The electric propagator is shown in Fig. > where the free massless propagator is also shown by the dashed 
line. Gluons at short distances, i.e., z < 6 — 1/T, have very similar behavior to the free propagator. However, at long 
distances, gluon propagators decrease more rapidly than the free one. This indicates that the electric screening mass 
does not vanish at all temperatures. 

Comparing data below and above Tc, we find that the dynamics of the gluon is completely different in the con- 
finement and deconfinement regions. The electric gluon mass in the confinement region becomes heavier; the electric 
gluon is completely screened in the confinement regions. This was first observed in Ref. Consequently, we 

cannot employ the assumption Eq. (|18() . On the contrary, the propagator in the deconfinement regions decreases 
exponentially even at long distances with a finite mass. 

At=0.05, a=1.0 



O o T/T =0.862 

a oTn-=i.05 

A A ^^=2.29 
TyT=4.97 
-- Free 



FIG. 6: The electric propagator in the confinement (circle) and deconfinement regions (other symbols) with Px(y) ~ '2TT/Nx(y)- 
All propagators are found to become massive compared with the free propagator (long dashed line). In the confinement regions, 
the propagator at long distances behaves like a very massive particle and vanishes, whereas the propagator beyond Tc has a 
finite mass. 



Similar behavior is seen in the magnetic parts in Fig. [7| although the long distance magnetic gluons have large 
errors here, all magnetic gluons are found to be massive at long distances in the confinement/deconfinement phase, 
and their effective mass in the confinement phase is heavier. We notice that the short distance behavior of magnetic 
gluons in the deconfinement regions looks unconventional. The magnetic gluon propagators follow a convex curve at 
short distances while the electric ones do not. The magnetic gluons behave as if they had an imaginary screening 
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mass or negative spectral function, which may be the reason why the magnetic mass is not screened at least in LOP 
calculation. 

As clearly seen in Fig. Hand as we discuss below, the effective mass of the magnetic mass is z dependent. We take 
the value around z = l/T as the electric case, since it is a relevant quantity at finite temperature screening. 



Ax=0.04, a=1.0 



z 

<N 0.01 



O-OT/T 


=0.867 


O OT/T 


= 1.15 




=2.69 


- - T/T 


=4.97 


-- Fl-ol 





{3 --i 



FIG. 7: The magnetic gluon propagator in the confinement (circle) and deconfinement regions (other symbols) with Pxiy) = 
27r/A^j;(y). Although the magnetic gluon propagators in this figure have large errors at long distances, they seem to have similar 
behavior to the electric part, except for the short distance behavior in the deconfinement regions. 




FIG. 8: Transverse gluon propagator, Grit) with p — {{2tt/Nx),0,0), on 48'^ x 64 at /3 = 6.8 (confinement phase). Solid line 
represents the free propagator |2l| . 



Gluons are essential ingredient of QCD but they are confined below Tc. As shown in Fig. |H1 their propagator is 
convex upward at several regions. This is possible only when the spectral function is not positive definite. This peculiar 
behavior was first observed in Ref. and confirmed in Ref. [33 ■ The feature does not contradict the fundamental 
postulate of quantum field theories, because gluons in the confinement region are not physically observable particles. 
Instead, this is a glimpse of the confinement mechanism in the infrared region which is still far from our understanding. 
Gribov's conjecture for the gluon propagator 

(19) 

vanishes at — 0, and its Fourier transformation to the coordinate space is not convex downward, |45| 

G(t) - -e''*'=°'''^cos(rtsin<?!) + 0), (20) 
r 

where r = (|p|'* + 6^)^/^ and (j) = 5 tan~^(&^/p*^). The upward convex shape of the gluon propagators in the decon- 
finement region may provide us with hints about the glue dynamics. 
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Below the critical temperature Tc, we obtain a similar result even for electric gluons at short distances. This 
seems natural since the perturbative argument in the confinement regions is generally not suitable and the confining 
correlation function would also give the negative spectral function. 

D. Mass as a pole 

To obtain the screening mass from the propagators, the following formula is used; 

Ge(™) ~ cosh[i;,(„)(p)(z - iV,/2)]. (21) 

We employ data for z ^ 1/T(= Nta)^ because the screening effect occurs at sufficiently long distances. All fittings 
are done from z = 6 to N^/l, and /NDF ~ 0(1), where NDF indicates the number of degrees of freedom. 

To obtain the final result at Ar — 0, we must extrapolate the data with respect to Langevin step width. The 
Runge-Kutta algorithm is applied to reduce the finite Langevin step Ar dependence js^l ■ We perform simulations for 
a set of parameters with At = 0.03 ~ 0.05. Table ITU and Fig. ^represent E{p) measured here versus Ar. The slight 
dependence of Ar enables us to use a linear function when fitting data. We finally obtain the mass from Eij)) by the 

T/T^=2.69 

1 I • , • , • , • 1 

o Magnetic 
0.8 ■ Electric 

° 0.6 ^ 

2 , . ■ 

CM 

3 ^— IE 

LU 

0.2 

I ' ' ' ' ' ' ' 

0.02 0.03 0.04 0.05 0.06 

Ax 

FIG. 9: Ar dependence of masses is slight. Finally to obtain a final value at Ar = 0, we use the linear function for the 
extrapolation. 

following lattice energy- momentum relation [l4j |4^ : 

sinh2:^=sinh^^+ysin^^. (22) 
2 2^2 ^ ^ 



E. Gauge invar iance 

The screening mass is physical and expected to be gauge invariant. However, since the gluon propagators defined 
by Eqs. H16|l and (|17|l are gauge dependent, it is important to check whether the screening masses obtained here are 
gauge invariant or not. In addition, since the magnetic mass cannot be defined by a perturbative calculation, it is 
particularly important to check its gauge dependence. In Fig. 1101 we show the gauge parameter a dependence of 
electric and magnetic masses. The gauge dependence of both screening masses is found to be very slight, namely, the 
result strongly suggests that they are gauge invariant and physical observables. 

F. Temperature dependence 

We study the temperature dependence of the screening mass in the range T/Tc — 1 — 6 which would be realized 
in high energy heavy-ion collision experiments such as RHIC or LHC 39j. Table UTTl and Fig. 1111 show electric and 
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TABLE II: This shows the typical example of the Langevin step dependence at T/Tc — 2.69. For all simulations for At — 0.03 — 
0.05, approximately 0.2 — 0.4 x 10® steps (measurements) are used after eliminating about 3000 — 5000 steps as thermalization. 
"p 7^ 0" means including the momentum Px(y) ~ '^n/N^f^y). We extrapolate these data to At = and then obtain mea(pfj_ = 
0) = 0.470(38) and mma{Pn = 0) = 0.199(52) using lattice energy-momentum relation. 



At 


Number of steps 


mma(p / 0) 


mea(p / 0) 


0.05 


275000 


0.400(08) 


0.544(22) 


0.045 


240000 


0.364(12) 


0.527(21) 


0.04 


380000 


0.381(11) 


0.508(21) 


0.035 


280000 


0.369(08) 


0.551(24) 


0.0325 


300000 


0.336(09) 


0.565(28) 


0.03 


320000 


0.367(14) 


0.551(24) 


0.00 




0.369(46) 


0.561(40) 



T/T =2.69, Ax=0.05, 



o Magnetic 
■ Electric 



O 0.6 



r OA - ^ ^ 

Q. ^ ^ 

LU 

0.2 



I ' ' ' ' ' ' ' 1 

0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 



FIG. 10: Gauge dependence for electric and magnetic screening masses. Gauge dependence of both screening masses is very 
slight. 



magnetic masses as a function of the temperature. The magnetic part definitely has nonzero mass in this temperature 
region. As T increases, both m^^^rnj/T decrease monotonically, and at almost all temperatures, the magnetic mass is 
less than the electric one, except very near Tc where the electric mass decreases very quickly as T approaches T^. 



G. Comparison with LOP and HTL resummation results 



We perform a fitting analysis for our numerical results using the following ansatz: 

^=Ce.9(T), !^ = a„.9^(r), (23) 

whose g dependence is predicted by the perturbative and 3D reduction analysis U 0| and we assume Cg and Cm are 
free parameters. In the following discussion, the data above T ~ 1.5Tc are used. Here we use the running couplings 
as 

^ ^'^^ 26olog(M/A) V 26o log(/./A) )' ^ ' 

and we set ^ = 2ttT, which is the Matsubara frequency as the renormalization point and A = 1.03Tc |^ as the QCD 
mass scale. 60 and bi are the first two universal coefficients of the renormalization group, 

60 = lliVc/48^^ 61 = (34/3)(7Vc/(167r2))2. (25) 
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TABLE III: Temperature dependence of the electric and magnetic masses which are extrapolated to the Langevin step Ar — 0. 



,/T T/Tc m,/T 



1.05 
1.15 
1.32 
1.36 
1.47 
1.74 
1.99 
2.29 



1.506(438) 
2.694(288) 
3.348(408) 
2.904(336) 
3.138(342) 
2.700(426) 
2.898(498) 
2.484(234) 



2.802(054) 
2.484(258) 
2.406(246) 
1.986(296) 
1.866(222) 
1.620(300) 
1.608(270) 
0.990(264) 



2.69 
2.99 
3.41 
3.88 
4.40 
4.97 
5.61 



2.820(228) 
2.892(234) 
2.190(450) 
2.292(222) 
2.598(168) 
2.310(084) 
2.106(390) 



1.194(312) 
1.590(318) 
0.960(168) 
0.852(318) 
1.134(414) 
0.804(638) 
0.486(336) 





Magncliu 
■ Ekclric 
















FIG. 11: Temperature dependence of electric and magnetic screening masses. The dotted line is fitted by the assumption, 
~ g^T. For the electric mass, the dashed and solid lines represent LOP and HTL resummation results, respectively. 



As a result, we obtain 



Cm 



1.63(3), xVNDF = 0.715, 
0.482(31), xVNDF = 0.979. 



(26) 



The scalings expected in Eq. (|23|l for electric and magnetic masses are found to work well. However, the magnitude 
of Ce is larger than Cf^^ = 1. On the other hand, for magnetic mass, a self-consistent inclusion technique in Ref. 
1^ gives rUg = 0.568g^T, which is close to our fitting result. 

The HTL resummation technique applying the free energy of hot gluon plasma has been widely discussed llC| . 
Rebhan gave a formula for the electric mass in the one-loop HTL perturbation theory and for the case of SU{3), 



'e,0 



1 



27r m 



e,0 



log 



2mp 



, viefl = gT. 



(27) 



Here we assume the magnetic mass to be of the order of g^. Substituting our fitted value for m™, we can solve the 
above equation iteratively. In Fig. ^] we show this HTL resummation together with the LOP result. The HTL result 
gives a better description than the naive perturbation, upon comparing with our numerical experiment. 

The electric mass was obtained also by a heavy qq potential from the SU{'i) Polyakov loop correlator at finite 
temperature in Refs. p^ll8|. Our results here are inconsistent with theirs, since the mass extraction from the heavy 
qq potential cannot be consistently performed due to ambiguity of its fitting assumption. In addition, a 3D reduction 



argument |12| has shown that me/ gT goes down when T increases, but even at T 
still about Zniefl. This observation agrees qualitatively with our analysis. 



lOOOA 



MS 



the electric mass is 
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TABLE IV: Data extracted from the small and large lattices in the same coupling regions, /3 — 7.5 and /3 — 8.0. 



13 At mma(p / 0) mea(p / 0) 


small lattice size 20'' 


X 32 X 6 


7.5 0.05 0.301(15) 
8.0 0.04 0.311(09) 


0.455(09) 
0.457(14) 


large lattice size 32^ 


X 48 X 6 


7.5 0.05 0.231(03) 
8.0 0.04 0.214(04) 


0.433(08) 
0.406(14) 



H. Higher T and finite size effect 

Although the main result in this paper is based on studies for the small lattice size 20^ x 32 x 6 as discussed in the 
previous section, we additionally perform the simulation on the large lattice 32^ x 48 x 6 to go to higher temperature 
regions and to check the finite-size effect of the screening masses. However, as the lattice size increases, the behavior 
of the long distance gluon cannot be controlled because of a large fluctuation. A typical result on the large lattice is 
shown in Fig. ^1 Even after the 0.3 — 0.4 x 10^ measurements, we could not determine the electric gluon propagator 
at long distances (z ^ 16), while the magnetic gluon is properly correlated. Nevertheless, provided that we adopt 
only the data for the intermediate regions above z = 6 until the disappearance of the propagator, we obtain similar 
results for the electric and magnetic masses as seen in Fis. 11 21 in the same temneratnre rerions. 

T/T=4.97, AT=0.05,a=1.0 



o G^(32'x48x6) 
O G (32^x48x6) 



FIG. 12: Typical electric and magnetic propagators on the large lattice 32^ x 48 x 6 (open), and the data for 20^ x 32 x 6 (band 
between dotted lines with error) are reproduced by using the values in Table ITTTl Both calculations give very similar value for 
electric and magnetic parts. 



Using the criterion describing above, we may consistently obtain both screening masses on the large lattice, and 
can argue that the magnetic mass has a stronger finite-size effect than the electric one. E{p) appearing in Eq.(17) are 
shown in Table HVI at T/Tc — 8.99 and 16.12. The momenta for the small and large lattices are Px = 2tt/Nx ~ 0.314 
and Px = 2'k/Nx ^ 0.196, respectively; namely, the result on the small lattice implies that the mag netic mass at high 
temperature {T/Tc > 5) seems to be going to zero, while on the large lattice it remains finite. [48j 

Although it is very difficult to measure the long distance gluon propagators, we can add the higher temperature 
results ( T/Tc = 8.99, 16.12 ) summarized in Table Ivl In Fig. E|we again fit the data including these new points. 
The fit for the large lattice data by Eq. (|23|) hence results in 



Cm 



1.69(4), 
0.549(16), 



X^/NDF = 
X^/NDF 



0.66, 
= 1.27 



(28) 



These results are shown in Fig. E| It should be noticed that Ce is the same value as given in Eq. (|26|l . while 
Cm ^ 10% larger on the large lattice and is very close to nig — 0.568(?^T calculated by the self-consistent inclusion 
technique in Ref. [ioj . 
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TABLE V: Simulation parameters and screening masses for the large lattice 32^ x 48 x 6. Lattice scales are estimated by Refs. 



P a-'[GeV] T[MeV] T/Tc m^/T m,/T 

7.0 7.64 1274 4.97 1.128(78) 2.556(156) 

7.5 13.8 2303 8.99 1.014(54) 2.178(144) 

8.0 24.7 4127 16.12 0.984(60) 2.256(120) 




T/T 

c 

FIG. 13: The temperature dependence including higher temperature points on the large lattice 32^ x 48 x 6. 



IV. CONCLUSIONS 

We have measured the gluon propagators and obtained the electric and magnetic masses by lattice QCD simulations 
in the quenched approximation for SU(3) between T = and BT^. Features of the QGP in this temperature region 
will be extensively studied theoretically and experimentally in the near future. 

Our screening mass studies are the first reliable measurement in SU (3) lattice calculation. We mainly investigate 
the temperature dependence for the electric and magnetic masses which do not vanish on 20^ x 32 x 6 lattices. In 
all temperature regions we find that the electric mass me is always larger than the magnetic one mm, except near 
critical temperature point. As the temperature goes down toward Tc, m^/T drops down quickly, while m,„/T is still 
going up. Consequently, using data above T/Tc ^ 1.5 we conclude that the scalings me ~ gT and m^ ~ g^T work 
well. Furthermore, a HTL resummation calculation has recently been developed and compared with nonperturbative 
lattice simulations. We have also compared our numerical results with LOP and HTL resummation and find a good 
improvement of the HTL electric mass. These comparison studies of SU{2>) screening masses qualitatively seem to 
agree with the case oiSU[2) 113. 

The electric masses obtained here are not consistent with those obtained by heavy qq potential calculations from 
an SU (3) Polyakov loop correlator at finite temperature in Refs. 0, 0| . In Ref. 13| , the authors did extensive 
analyses with three different temporal extents and two different gauge actions, obtaining a very reliable potential as 
a function of the temperature. They observe that the potential above Tc cannot be described properly by the leading 
order perturbation calculation up to a few T^. They exclude the two-gluon exchange as the dominant screening 
mechanism, and suggest that some kind of one-gluon exchange may describe the potential effectively as a result of 
the complex interaction, and that at about (1.5 — 3)rc a mixture of one- and two-gluon exchange may explain the 
behavior. Therefore, due to the ambiguity of the fitting assumptions, it is not clear whether we can compare our 
screening masses directly with those obtained by the potential calculation. 

In order to investigate the nature of the QGP, especially the excitation modes in the plasma, Datta and Gupta 
recently calculated glueball masses at finite temperature and made an interesting observation. They measured the 
screening masses of A^^ (scalar) and (glueball), which allow two- and three-gluon exchange, and their ratio ~ 1.7 
is near 3/2. The mass is twice that obtained by Kaczmarek et al, and shows similar temperature dependence. 
There are now several nonperturbative methods to study QGP: our direct measurement of the gluon propagators. 
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glueball screening masses, and Polyakov line correlators. These analyses strongly suggest that the QGP above Tc is 
far from a free gas and has a nontrivial structure. Much more detailed analyses in future are highly desirable. 

The screening mass on the lattice is extracted from the gauge dependent propagator, and the magnetic mass is 
not well defined in perturbation theory. We have nonperturbatively confirmed the gauge invariance of both screening 
masses. In Ref. 20] it was reported that the SU{2) magnetic propagator exhibits a complicated gauge dependent 
structure at low momentum. Therefore, since the gauge dependence for the screening masses is investigated within 
Lorentz-type gauge fixing based on stochastic gauge quantization in this study, we plan to extend our analysis to a 
simulation with Coulomb-type gauge fixing. 

We have seen a qualitative difference of the gluon dynamics between the confinement and deconfinement phase by 
direct propagator measurement. The electric and magnetic gluons in the confinement phase indicate a very massive 
particle behavior, while after the phase transition, they have a finite mass. In addition, in the deconfinement phase, 
the magnetic gluon at short distances seems to be still in the confinement phase. This may be related to the difficulty 
of the perturbative argument for spatial gluon components and the fact that a magnetic Wilson loop gives non-zero 
spatial string tension even at high temperatures |43| . 

The magnetic mass has been the subject of many discussions. Perturbatively, it is difficult to handle. To our 
knowledge, there is no complete perturbative calculation which is free from any assumption or model. The naive 
expectation is that it vanishes, but it is necessary to have a finite value as a cutoff factor in the infrared regime. On 
the other hand, the finite spatial string tension even at T > Tc indicates "confinement" in the magnetic sector. 

Our screening gluon magnetic propagators indeed show nontrivial behavior. At very short distances, this may be 
consistent with massless behavior, but at finite distance we cannot fit them by a simple ansatz. The behavior is 
distance dependent. Confinement is a long range property, and the propagators there drop. Therefore it will be a 
very interesting task in future to investigate the magnetic propagators at very long distance. 

We calculated the gluon propagators on the larger lattice 32^ x 48 x 6 and observed a large error and strange 
behavior at long distances. Nevertheless, the screening masses were estimated, and we find that the magnetic mass 
is sensitive to the lattice size effect. Thus on a too small lattice we cannot deal with the magnetic mass consistently. 
Moreover, the simulations even at higher temperature T/Tc 9 and 16 show that nonperturbative results are far 
from LOP. This observation is compatible with that of Ref. [l^ IT^ . 

For quantization with gauge fixing, stochastic gauge fixing was adopted. We think the stochastic gauge fixing has 
better features to reduce some of the difficulties of nonperturbative gauge fixing. It is consequently possible to do 
a practical simulation of gluon screenings effectively. However, we also see that the gluon propagators have large 
fluctuations and unexpected behavior at long distances and may need further calculations. 

The color screening data we obtained here are useful information for QGP phenomenology, for instance, jet quench- 
ing or the heavy quark potential. We plan to study the qq as well as the qq potential relating a baryon bound state, 
a nonperturbative QCD vertex calculation, quark propagators, etc., using stochastic gauge fixing, which will help us 
to understand QGP. 
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